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ABSTRACT 

We present sensitive near-infrared spectroscopic observations for a sample of five z ~ 6 quasars. 
These quasars are amongst the most distant, currently known quasars in the universe. The spectra 
have been obtained using ISAAC at the VLT and include the C IV, Mg II and Fe II lines, which are 
powerful probes of the chemical enrichment and the black hole masses in these objects. We measure 
the Fe II/Mg II line ratio, as an observational proxy for the Fe/a element ratio. These elements 
are produced by different types of supernovae and their abundance ratio can therefore serve as a 
cosmological clock. We derive a ratio of 2.7±0.8 for our sample, which is similar to that found for 
lower redshift quasars, i.e., we provide additional evidence for the lack of evolution in the Fe II/Mg II 
line ratio of quasars up to the highest redshifts. This result demonstrates that the sample quasars must 
have undergone a major episode of iron enrichment in less than one Gyr and star formation must have 
commenced at z > 8. The linewidths of the Mg II and C IV lines give two estimates for the black hole 
masses. A third estimate is given by assuming that the quasars emit at their Eddington luminosity. 
The derived masses using these three methods agree well, implying that the quasars are not likely to 
be strongly lensed. We derive central black hole masses of 0.3 — 5.2 x 10^ M©. These include the lowest 
black hole mass ever measured at z 6, suggesting that we probe a more typical quasar population 
(with lower masses/luminosities) than examined before. We use the difference between the redshift of 
Mg II (a proxy for the systemic redshift of the quasar) and the onset of the Gunn Peterson trough to 
derive the extent of the ionized Stromgren spheres around our target quasars. The derived physical 
radii are about five Mpc. Using a simple ionization model, the emission of the central quasars would 
need of order 10® — 10* yr to create these cavities in a surrounding intergalactic medium (IGM) with 
a neutral fraction between 0.1 and 1.0. As the e-folding time scale for the central accreting black hole 
is on the order of a few times 10^ years, it can grow by one e-folding or less within this time span. 
Subject headings: Galaxies: high-redshift — Galaxies: fundamental parameters — Galax- 
ies: formation — Galaxies: evolution — Quasars: individual: J000552.34- 
000655.8, J083643.85+005453.3, J103027.10+052455.0, J130608.26+035626.3, 
J141111. 29+121737.4 



^ Based on observations carried out at the European South- 
ern Observatory, Paranal, Chile under program Nos. 267.A-5689, 
069.B-0289, 071.B-0525, 074.A-0447 and 076.A-0304. 

^ German Science Foundation's Collaborative Research Center 
439 fellow. 



1. INTRODUCTION 

Since the advent of the Sloan Digital Sky Survey 
fSDSS. lYork et al.l l20GG) and other large muhi-band sur- 
veys, more than twenty quasars at z > 5.7 hav e been 
discovered (mostly by the SDSS: iFan et al.l[2000l l200lL 
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20031 12004 I2006d and three in other surveys : iGotol 



2006HMcGreer et al.l[200llVenemans et al.ll2007f ). These 



quasars provide the best probes of the early growth of 
supermassive black holes and serve as signposts for over- 
densities in the universe only one Gyr after the Big Bang. 
High redshift QSO studies also provide imp ortant in- 
sight into the reionization (|Fan et al.ll2006af i of the in- 
tergalactic medium (IGM) that took place in the early 
Universe and into the relation between t he formation of 
early galaxies and black holes (|Fanll200 6[). 

Studies of the full SPSS quasar sample from = 
to z = 6 (jRichards et al.ll2002at [Schneider et"all l2005( l 
show a strong evolution in number density of the quasar 
phenomenon with redshift (the comoving density of lu- 
minous quasars a t z ^ 6 is about forty times smaller 
than at z ~ 3, [Richards et al.l l2006ct) . but there is 
no evolution detected in th e intrinsic spectral proper- 
ties of these quasar phases ()Fanl l2006l ). The high den- 
sity peaks in the dark matter distribution, which host 
QSOs, are ex pected to be much rarer in the early uni- 
verse (iShen et a l. 2007), while the immediate environ- 
ment of quasars matures apparently very early on. Blue- 
ward of the Lya line, the observed spectra of SDSS 
quasars do change strongly as a function of redshift, as 
the absorption by intervening neutral hydrogen increases 
with redshift up to the point where the continuum emis- 
sion is completel y absorbed (complete G unn-Peterson 
(GP) absorption, iGunn fc PetersonI 119651) . The detec- 
tion o f deep GP troughs in 23 z ~ 6 QSOs (jFan et al.l 
l2006cf ) indicates that the IGM is significantly neutral 
(riHi/'T-H ^ 10~^ — 10~^) at these redshifts. However, the 
large dispersion of IGM properties measured along differ- 
ent lines of sight strongly suggests that the reionization 
process is more complex than a quick phase tr ansition 
over a narrow redshift range C Fan et al.l[2006bh . This 
finding is supported by the results of the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP. iBennett et al.]l2003f l 
mission, which suggest that the universe may have been 
~ 50% neutral at redsh ift z > 10, but r eionization did 
not start before z ~ 14 ([Page et al.ll2007t l. QSOs acting 
alone cannot have reionized the hydrogen in the IGM 
prior to z ~ 4 fe.g. . 'Meiksinl l2065l ) . leaving star-forming 
galaxies as the most likely candidates to ionize the IGM 
at z > 6. 

The strong emission lines excited by the nuclear source 
of a galaxy can be used to infer various physical and dy- 
namical properties of the circumnuclear gas and therefore 
its power source, the nuclear black hole. For QSOs at z > 
5.7, the strong emission lines redward of Lya are shifted 
into the near infrared. Several near-infrared (NIR) spec- 
troscopy studies of QSO sam ples including objects at z ^ 
6 have been carr icd out (e.g.,[Ma iolino et al 2001, '2003; 
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studies show that the Fe II/Mg II line ratio appears to 
be similar at all redshifts, up to and including z ^ 6 
(some studies even report higher values than found in 
the local universe). Assuming that this ratio is directly 
related to the Fe/a abundance, this suggests that all ob- 
served quasars have undergone a major episode of iron 
enrichment early on, at least in their nuclear region. We 
do note one caveat here: at high redshifts, we do not 
know about the metallicity of less luminous QSOs as, so 



far, only the brightest objects have been studied. 

The Fe to Mg abundance ratio, and its observational 
proxy, the Fe II/Mg II line ratio, can be considered a 
cosmological clock. Both elements are produced in su- 
pernova explosions, but while Fe is produced by Type la 
supernovae (SNe), which have relatively low mass pro- 
genitors (white dwarfs in binary systems), Mg is pro- 
duced by Type II SNe, which have high mass progen- 
itors. Mg therefore appears almost instanteneously af- 
ter initial star formation while the Fe production starts 
only later. The ratio of Fe to Mg is predicted to build 
up quickly in the first 1 to 3 Gyr and then level off to 
the v alue presently obse rved in the solar neighbourhood 
(e.g., lYoshii et al.lfl99l . More recent studies, however, 
show that this may be true for our Milky Way and sim- 
ilar galax ies, but can be as short a s 0.3 Gyr for elliptical 
galax ies ()Friaca &: TerlevichI 119981 : iMatteucci fc Recchil 
l2001h . as the time of maximum enrichment is a strong 
function of the adopted stellar life times, initial mass 
function, and star formation rate. 

Evidence for enrichment on kiloparsec scales is also 
available from the detection of spatially resolved carbon, 
oxygen and du st in the int erstellar me dium around a 
QSO at z = 6 4 (IWalter et al . 2003; Bcrtol di et al.ll2003al : 
IWaltcr et al.' 2004; Maiohno et al.ii2005i) . ISimcod pOOef ) 
and Ryan- Weber et al. (200f|) have observed two z 6 
QSOs in the NIR at moderately high spectral resolution 
and have found several intervening C IV absorption sys- 
tems, suggesting that f^civ does not evolve with redshift. 
This therefore indicates that even at Mpc scales a large 
fraction of intcrgalactic metals may already have been in 
place at z > 6. 

As the width of the observed emission lines is believed 
to represent (at least in part) ordered motion of gas in 
the very center of the quasars, these can be employed to 
estima te the mass which drives this motion (see, for ex- 
ampl e, iMcLure fc Dunlopll2004l : IVestergaard fc PetersonI 
l2006f ). The width of the Mg II line is used by several 
authors to derive masses in the range 1 — 6 x 10^ M0 for 
the central black hole in distant QSOs. For the redshift 
record holder, quasar SDSS J1148+5251 at z = 6.42, this 
indeed indicates that the quasar is accreting at the maxi- 
mal allowable rate for a black hole, ado pting the Edding- 
ton limit criterion (jWillott et al.l[2003l ). This quasar has 
also been detected at mm wavelengths, provi ding impor- 
tant in formation about its dust content (Be rtoldi et al.l 
[20033), the. mass of the Ho gas reservoir in the host 
galaxy (I Walter et al.ll2003D. the densitv and temperature 
of CO gas (iBertoldi et al.ll2003H ) and its dynamical mass 
() Walter et al.ll2004l ). The estimated dynamical mass of 
~ 5 X 10^° M0 is inconsistent with a ^ 10^^ M© stellar 
bulge predicted if the Mbh ~ fbuigc relation (which seems 
to be valid at lower redsh ifts, e.g.. iGebhardt et al.ll2000l : 
iFerrarese fc Merrittl l2000t ) were to hold at z ~ 6. This 
indicates that black holes may form pri or to the assem- 
bly of stellar bulges ([Walter et al.ll2004D . The detection 
ofJC II 158/im] emission in this QSO bv ^Maiolino et all 
()2005| ) shows that i ts host is un d ergoin g an intense burst 
of star formation. IWang et"al] (|2007l ) present a sample 
of eighteen z '--^ 6 QSOs observed at submm wavelengths. 
Eight of these are detected at mJy sensitivities, strength- 
ening the idea that massive starbursts exist in quasar 
host galaxies. 

Some of the above mentioned NIR studies have been 
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hampered by low S/N and/or limited wavelength cover- 
age of the obtained spectra, or included only very few 
z ~ 6 QSOs. With the goal of measuring the physi- 
cal properties from high-quality observations of a larger 
number of z ~ 6 quasars in a consistent way, we have ob- 
served a sample including all published QSOs at redshift 
z > 5.8 and below declination S = +15°, i.e., those acces- 
sible from the four 8m apertures that constitute the Very 
Large Telescope (VLT). Our NIR ISAAC spectroscopy in 
several bands from 1.0 to 2.5 /im covers the C IV, Mg II 
and Fe II lines. The sample includes objects with lower 
luminosities than observed in earlier studies and thereby 
helps to constrain the properties of presumably more typ- 
ical QSOs at these redshifts. 

This paper is divided in the following parts: in Sec. [2] 
we describe the sample, the observations and the data 
reduction we have carried out. In Sec. [31 the results 
are shown, including an explanation of the fitting pro- 
cedure, followed by the measurement of the properties 
of the spectral features. In Sec. IH the interpretation 
of these results is presented, amongst others, in terms of 
metal enrichment and black hole mass. Finally, we briefly 
summarize our results in Sec. [5l We assume the follow- 
ing ACDM cosmology throughout the paper: Hq = 70 
km s-i Mpc-i, = 0.3, Qi, = 0.04 and f^A = 0.7. 
These parameters are al so use d in papers written by 
iMcLure fc DunlopI (|2004f l and IVestergaard fc Petersonl 
(I2006D . from which we use several relations. 

2. SAMPLE DEFINITION, OBSERVATIONS AND 
DATA REDUCTION 

Of the eleven quasars at z > 5.8 known in 2003, we 
have selected those that are easily observable from the 
Southern Hemisphere, i.e., at declinations lower than 
15°. This resulted in a sample of five targets with mag- 
nitudes 18.7 < z^g < 20.5 and redshifts from discovery 
papers in the range 5.8 < z < 6.3 (see Table[l]for names, 
redshifts, magnitudes and discovery papers). 

The obs ervations were carried o ut with the ISAAC in- 
strument (jMoorwood et al.lll998D on Antu (VLT-UTl) 
in low resolution (LR) mode, using the 1024x1024 
Hawaii Rockwell array of the Short Wavelength arm. 
ISAAC'S LR grism covers four bands, called SZ, J, 
SH, and SK (referred to as Z, J, H and K in this 
paper). Depending on the order selection filter chosen, 
this results in spectral resolutions of 550, 500, 500, and 
450, respectively, in combination with the 1" slit that 
we have employed (for the Z- and J-band observations 
of J1030-f0524 and J1306+0356, respectively, we used 
a narrower 0'.'6 slit resulting in a slightly higher resolu- 
tion). For each QSO, the C IV and Mg II lines were 
observed. Depending on the redshift, these lines are in 
the Z- or J-band (for J1030-)-0524), and X-band, re- 
spectively. For J0836-I-0054, a spectrum in the iJ-band 
was also obtained. Table [T] summarizes the bands chosen 
and the exposure times for each objects. The observa- 
tions were carried out in service mode and were collected 
over several years: some of the data were taken as early 
as May 2001, while the bulk was obtained between Octo- 
ber 2004 and March 2006. Exposure times were typically 
2.0 hours in Z-band and 6.0 hours in ii'-band, but vary 
from object to object (see Table [1]). 

In most cases, twenty frames of 120 seconds were ex- 
posed in one hour (or observation block), in ABBA se- 



quences with large 20" to 30" offsets (dithering) and 
small random offsets (jittering to avoid bad pixels, within 
a box of 4" to 18"). A telluric standard was observed 
each night (before or after a quasar was observed) at an 
airmass typically within 0.2 of the airmass of the QSO 
observations. 

The data were reduced using ESO's pipeline software 
and IRAF. Although ESO delivers the observations to- 
gether with pipeline reduced products, we have rereduced 
the original frames using the newest version of ESO's 
pipeline software (for ISAAC, version 5.4.2) to ensure 
that all spectra were reduced in the same way. In partic- 
ular, care was taken that the pipeline reduction was car- 
ried out using sky emission lines to determine the wave- 
length calibration (this was not always the case in the 
ESO products). Given the short integration times of the 
telluric stars in the Z-band, the imprints of the sky lines 
were not strong enough to be used as a wavelength indi- 
cator, and we thus used the wavelength calibration given 
by arc lamps observed the following day. The zeropoints 
of the wavelength scales of the Z-band spectra were then 
corrected by measuring the position of a single bright sky 
line and, in a few cases, of the hydrogen Paschen lines 
in the continua of stars serendipitously observed in the 
same slit as the target. 

Rerunning the pipeline produces wavelength calibrated 
co-added images of the individual frames observed in 
one observation block. Subsequent reduction was car- 
ried out within IRAF. Using IRAF's apall task, one- 
dimensional spectra were extracted from these images. 
The extraction width was either 5 pixels (0'.'75) or the 
seeing FWHM. The spectra were then corrected for tel- 
luric extinction using IRAF's telluric task (by dividing 
the observed QSO spectrum by the observed star spec- 
trum), after removal of instrinsic spectral absorption fea- 
tures from the observed star. To get the correct slope and 
flux calibration, the spectra were multiplied by an artifi- 
cial black body spectrum of the standard star based on 
its spectral type and magnitude (in the observed band if 
available from the literature, but otherwise in the clos- 
est band for which a magnitude was available). Then, 
the individual one-dimensional spectra were averaged to 
form a single spectrum per quasar per band. The flux 
calibration of these spectra was checked by comparison 
with broad band magnitudes. As most of the necessary 
broad band magnitudes were not available from obser- 
vations, we fitted the type 1 quasa r template SED con- 
structed by [Richards et al.l (l2006bD to the avail able NIR 
and MIR photometry from iJiang et al] (|2006f) and de- 
rived Z, J and K flux densities from the resulting nor- 
malized SEDs. As the final step in the data reduction, we 
corrected the absolute fiux calibration by factors within 
the range 0.75 — 2.2 to match these derived broad band 
fluxes. 

3. RESULTS 

3.1. Composite spectrum 

For illustrative purposes, we have constructed a com- 
posite spectrum by shifting the five QSO spectra to their 
restframe wavelength and averaging them, thereby ob- 
taining a spectrum covering a larger wavelength region 
than any of the single spectra (Fig.[T]). The spectra were 
normalized to all have the same mean fiux density in the 
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TABLE 1 
VLT ISAAC Observations 



FuUname^' Rcf^ z z^'g t^'P^'^ t"/P''* t^P'<* ^^^'"^ 

[h] [h] \\ \\ 



SDSS J083643.85 


+ 005453.3 


(2) 


5.82 


18.74 


1.7 




3.2 


4.0 


SDSS J000552.34 


- 000655.8 


(3) 


5.85 


20.54 


3.3 






5.3 


SDSS J141111.29 


+ 121737.4 


(3) 


5.95 


19.64 


2.0 






6.0 


SDSS J130608.26 


+ 035626.3 


(2) 


5.99 


19.47 


1.6 






4.6 


SDSS J103027.10 


+ 052455.0 


(2) 


6.28 


20.05 




1.6 




12.0 



^ Full names arc listed here. In the main text and following tables, only appropriate ab- 
breviations are used, e.g., J0836+0054 for SDSS J083643.85 + 005453.3.^ Discove ry papers, 
from which redshifts and magnitudes presented in the table have been taken: {V\ iFan et al.l 
120001 . (2) IFan et"ani2001| . (3) IFan et"ani2004i '^ Magnitudes in ^'-band from discovery papers 
in asinh AB system. Magnitude errors are between 0.05 and 0.10.*^ Exposures time in hours 
for the bands indicated. 



restframe wavelength range 3000 - 3200 A. As absolute 
flux reference in this region was taken the spectrum of 
J0836+0054. Several emission lines are visible, most no- 
table Lya and Si IV/0 IV] in the observed-frame optical 
domain, and C IV, Mg II and Fe II in the observed-frame 
NIR domain. To measure the properties of these lines, 
we must determine the underlying continuum, and, for 
Mg II, we also need to determine the Balmer pseudo- 
continuum and the Fe II emission line forest. This spec- 
tral decomposition is described in the following section. 

3.2. Decomposition of emission lines and continuum 
3.2.1. Continuum 

The dominant component in the spectra is the con- 
tinuum, which is usually modelled as a power law {F\ cx 
A"). For the quasar composite spectrum obtained from a 
sample of over 2200 SDSS spectra within a redshift range 
0.044 < z < 4.789 (see Fi g. H), the slo pe was found to 
be a\ = —1.6 (jVanden Berk et al] 1200 ll). For individual 
quasars this value varies, e.g., [Richards et aH (|2003) find 
values between -1.75 and -1.24 for most of the 4576 SPSS 
quasa rs at 0.3 < z < 2.2 in their sample. [Pentericci et al.l 
(|2003') find a value of -1.43 with a scatte r of 0.33 for 45 
SPSS q uasars at 3.6 < z < 5.0, while iMaiolino et al.l 
(|2004aD find power law slopes between -1.35 and -2.10 
for a sample of eight quasars at 4.9 < z < 6.4. For spec- 
tra with short wavelength coverage and many emission 
lines (e.g., Mg II, Fe II), the determination of the power 
law slope is made difficult by the fact that almost no 
part of the quasar spectrum is not also covered by one 
or several (blended) emission lines. We, therefore, do fit 
a power law to be able to measure the flux in the emis- 
sion lines superposed, but note that the measured slope 
has a large error and might not be representative for the 
overall continuum shape of the quasar. 

3.2.2. Balmer pseudo-continuum 

A second component which has to be determined and 
subtracted in order to make proper fits to the Mg II 
and Fe II emission lines, is the forest of hydrogen emis- 
sion lines, which form the so-called Balmer pseudo- 
continuum. 

To model the Ba lmer continuum, we have followed 
iDietrich et all (|2003f ). who assume gas clouds of uniform 
temperature (Tg = 15000 K) that are partially opti- 
cally thick. The Balmer continuum spectrum below the 



Balmer edge (Abe = 3646 A) can be described by 

^^Bac ^ F^^Bx{Te){l - e-^BE(A/ABE)=')^ X < Abe (1) 

where B\{Te) is the Planck function at the electron tem- 
perature Te, tbe is the optical depth at Abe, and F^^ is 
the n ormalized flux density at the Balmer edge (jGrandil 
Il982f) . In general, the normalization is estimated at 
A ~ 3675 A since no Fe II emission is present at that 
wavelength. However, in our spectra, this region is cither 
absent or has a low signal-to-noise ratio. Therefore, the 
normalization is fixed to the strength of the power-law 
continuum during the fit and the optical depth is fixed 
at Tbe = 1 (we found minimal differences in the Balmer 
continuum when we varyied tbe between 0.1 and 2.0). 

3.2.3. Fe II template 

The Fe II ion emits a forest of lines, many of which are 
blended. These are usually fit jointly by using an iron 
emission line template. Several templates are available in 
the literature, based on high signal -to-noise ratio spectra 
of ind ividual active galaxies (e.g.. iVestergaard fc Wilkei 
1200 ID . on com posites of many quas ar spectra (e.g., from 
the LBQS bv iFrancis et~al] [l99ll and from the SPSS 
by iVanden Berk et aL 2001), or on theore tical models 
(e.g.. iVerner et al.l |1999.: Sigut fc Pradha^ [200l . Be- 
fore fitting our data, we have compared two templates 
applicable to high redshift quasar spectra: one derived 
from the nearby narrow line Seyfert 1 g alaxy , I Zwicky 
1, pubfished by IVestergaard fc Wilked ([2001') and one 
we derived from the com posite SPSS QSO spectrum 
(jVanden Berk et al.ll200l[) . 

/ Zw 1 template by Westeraaard & Wilkei 

IVestergaard fc Wilked have constructed an Fe II and 
Fe III template based on Hubble Space Telescope spectra 
of the narrow line Seyfert 1 galaxy, I Zwicky 1 (I Zw 
1, PG 0050-^124, z = 0.061). Its narrow intrinsic fines 
(~ 900 km s~^) and its rich iron spectrum make the tem- 
plate particularly suitable for fitt ing the spectra of ac- 
tive galactic nuclei and quasars. IVestergaard fc Wilkes! 
fitted and subtracted a power law continuum, and ab- 
sorption and emission features from other elements than 
iron from the I Zw 1 spectrum. An Fe III emission model 
was subtracted from the residual to create a pure Fe II 
template. In the process, the Mg II AA2795,2803 emis- 
sion was fit with two double Gaussian components 400 
km s^^ apart. After subtraction of this Mg II fit, no 
iron emission is left in the 2770 — 2820 A range, although 
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Fig. 1. — Composite spectrum of the five 2; ~ 6 QSOs in our sample at rest frame wavelength. Also displayed is a histogram showing the 
number of spectra contributing at each wavelength. Strong emission lines are identified. The depression in flux below the hya line shows 
where the GP trough is located (although this is much more pronounced in individual spectra). 



from theoretical considerations, it should not be entirely 
absent. Thus, subtraction of their Fe template would 
leave some unidentified iron emission around 2800 A, 
leading to an overestimate of the flux of the Mg II line. 
We have therefore add ed a constant flux density to the 
IVestergaard fc Wilkes! template of ~20% of the mean flux 
density of the template between 2930 and 2970 A, as il- 
lustrated in Fig. [3] (rough ly consistent with the value 
envisaged by Vcstcrgaard fc Wilkesll200ll) . The justifica- 
tion for this operation co mes from the theor e tical Fe II 
emission line strength by iSigut fc PradhanI (|2003l see 
their Fig. 13). 

SDSS QSO composite template by I Vanden Berk et al\ 

From the composite SDSS spectrum, smoothed to 
a resolution of 7 A, comparable to the observed rest 
frame resolution, we have sub tracted a power law contin- 
uum with the slope given bv lVanden Berk et all (f200l[ ). 
We have normalized the power law to the measured 
values in the 4200 - 4230 A region (as specified in 
IVanden Berk et ani2001h . In addition, a Balmer contin- 
uum was created and subtracted, normalized to the mea- 
sured value at 3675 A and with optical depth tbe = 1.0 
(again, only minimal differences were found when chang- 
ing Tbe within reasonable limits). The result is shown as 
an inset in Fig. O 



Comparison of the I Zw 1 and SDSS QSO composite 
templates 

We have subsequently compared the 

IVestergaard fc Wilkes! and the SDSS templates (see 
Fig. [3]). The figure shows that the two templates are 
very similar, apart from the Mg II line which we did not 
try to subtract in the SDSS template and an Fe III UV47 
line at A = 2418 A. We have the refore decided to 
use the IVestergaard fc Wilkes! (|2001f ) template for the 
analysis that follows, as it allows us to fit the Mg II line 
simultaneously. 

3.3. Fitting procedure 

Fitting was carried out with the MPFIT'^ package which 
performs the Levenberg-Marquardt least -squares mini- 
mization. Nominal errors of the fitted parameters were 
computed by multiplication of the formal one-sigma 
errors computed from the covariance matrix with the 
square root of the reduced x^, which was between 1.0 
and 1.5 in all cases. 

The if -band spectra were fitted in two steps (see 
Fig. 13]). First, the part of the spectra corresponding to 

^ M PFIT is an IDL routine written a nd distributed by C. Mark- 
wardt ) |http: //cow .physics .wise ■ edu/$ \sim$craigm ). 
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Fig. 2. — The SPS S median quasar composite spectrum 
llVanden Berk et al.l 1200 11). shown with an arbitrary flux scaling. 
The dashed line is a power law contiimum {Fx = X") with a slope 
a = —1.54 given by IVanden Berk et al.. and a normalization fit 
by us. The inset shows a zoom in the region of the Mgll and 
part of the Fe II emission. Also shown is the power law continuum 
and the co— added power law and Balmer continua, which are sub- 
tracted from the smoothed composite spectrum to form the final 
Fell emission template, shown in grey [blue]. [See the electronic 
edition of the Journal for a color version of this figure.] 
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Fig. 3.— Our modified Fe II |Vestergaard & Wilkefj lf20M > tem- 
plate (solid line, smoothed here for display purposes), compared 
with our template derived from the S DSS composite QSO spec - 
trum (indicated by the dashed [red] line. IVanden Berk et al.]|2001l) . 
The templates differ significantly only around 2418 A, where the 
Fe III UV47 line is located and around 2800 A, where the promi- 
nent Mgll line is located. These lines are not present in the 
IVe stergaard & Wilkes template which contains only Fe II emission. 
[Sec the electronic edition of the Journal for a color version of this 
figure.] 

rest-frame wavelengths in the range 2848 to 3081 A was 
fit simuhaneously by the power-law continuum, Balmer 
pseudo-continuum and iron template. Here, the Balmer 
continuum normalization was fixed to the power-law con- 
tinuum strength at 3600 A, which leaves only three free 
parameters being fit at this stage: the slope and nor- 
malization of the power-law continuum, and the normal- 
ization of the iron template. The iron normalization is 
mostly determined by the emission feature at ~ 2950 A. 
We did not attempt to fit the redshift of the iron tem- 
plate. Instead it was kept fixed (to the Mg II redshift) 



and only the template normalization was varied (a few 
iterations were necessary to obtain this redshift). The fit 
was then extrapolated to the full spectral range of the 
i^T-band spectrum and subtracted from the data. The 
residual still contains the Mg II line which was subse- 
quently fit with a Gaussian function, having three free 
parameters: central wavelength, width and normaliza- 
tion. The Fe II flux was computed by multiplying the 
derived iron template normalization with the total Fe II 
flux derived from the iron template by integrating it over 
the wavelength range 2200 < A < 3090 A. Here, the as- 
sumption is made that the model (the iron template) 
represents the observed data well. The obtained dur- 
ing the fltting procedure gives an indication for this rep- 
resentativeness, but only for the relatively small range 
2848 < A < 3081 A. IVestergaard fc Wilkes! (|2001l ) state 
that the iron template derived from I Zw 1 is generally 
applicable to high redshift quasar spectra. Our compar- 
ison to the SDSS template confirms this. To account for 
this uncertainty, however, we multiply the nominal er- 
ror in the Fe II flux by a factor two to obtain the error 
displayed in Table [H 

The Z and J-band spectra were fitted in one step (see 
Fig.O. The C IV line was fit with a Lorentzian function, 
while the underlying continuum was fit with a power law. 
The function fit has therefore five free parameters: the 
power law slope and normalization, and the Lorentzian 
function's normalization, width and central wavelength. 
J1411-I-1217 shows associated absorption in the C IV line, 
which was fit with a Voigt profile, adding four more free 
parameters. The fitting parameters are summarized in 
Table [3 

4. DISCUSSION 
4.1. Metal enrichment 

As described in Sec. [l] a measurement of the Fe to 
Mg ratio at z ~ 6 may serve as an approximate indica- 
tion of the onset of star formation in the highest redshift 
quasars. The measured Fe II/Mg II line ratios of our five 
z 6 QSOs are shown in Table [2] The ratios lie within 
the range 2 — 5. 

We have compared the Fe II/Mg II line ratio measured 
in our study to those of other studies at similar and lower 
redshifts (summarized in Fig. [6]). It is, however, difficult 
to make a direct comparison: different studies use differ- 
ent wavelength regions to estimate the underlying con- 
tinuum and Fe II feat ures and use different r nethods to fit 
the spectral features. Ilwamuro et al.l (|2002| ) try to over- 
come this problem by applying the same fitting algorithm 
to new spectra at 4.4 < z < 5.3 and archival spectra at 
lower redshifts down to z = 0.03. Although the scatter 
within their redshift bins is large, they measure a median 
value of Fe II/Mg II at z ^ 5 which is about 50% higher 
than at z ~ 1.5. Several other authors have studied 
the Fe II/Mg II fine ratio for QSOs at higher redshifts: 
5.7 < z < 6.3 (iFreudling et al.| [200l. 3.0 < z < 6.4 
(iMaiolino et al.ll2003f ) and 6.1 < z < 6.3 (jlwamuro et al.l 
|2004). The latter paper summarizes these data and con- 
cludes that the upper envelope of the Fe II/Mg II dis- 
tribution decreases towards higher redshift (z > 3) while 
the scatter among individual objects increases. This may 
indicate that some objects are observed such a short time 
after the initial starburst that the Broad Line Region 
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Fig. 4. — iC-band spectra of the five 2 ~ 6 QSOs, covering the wavelength region that includes Mgll, shown here smoothed over 5 pixels 
(or 0.0036 ^m). The observed wavelength is indicated below the spectra (in microns), while the rest frame wavelength is indicated above 
the spectra (in A). Several features of the fit are indicated on top of the spectra (dark [black] solid line): the power law continuum (grey 
[red] solid line), the Balmer pseudo-continuum (grey [red] dashed line), Vestergaard Fell template (light [green] solid line, thick where it 
is actually fitted to the spectrum), the Gaussian fit to the Mgll line (dark [blue] thick line). The dashed grey [blue] line shows the noise 
level per pixel (smoothing is not taken into account here). [See the electronic edition of the Journal for a color version of this figure.] 



(BLR) is not yet fully enriched with iron. Fig. [S] shows 
that our measurements confirm the lack of evolution in 
the Fe Il/lVIg II line ratio observed by other authors, in- 
dicating early enrichment of the BLR. As the age of the 
Universe at z = 6.0 is only 0.9 Gyr, an enrichment time 
of 1 Gyr needed to reach the observed ratio according to 
some models (e.g.. lYoshii et al.lll998l l seems to be ruled 
out. Assuming an enrichme nt time of 0.3 Gyr appropr i- 
ate for an elliptical galaxy flVl atteucci fc Recchil l200ll ). 
t he onset of s t ar for mation would have been at z > 8. 
IWills et all (|l985l) discuss how well the Fe U/Mg II 



line ratio represents the Fe/Mg abundance ratio, as the 
region where Fe II is produced is more extended than 
the Ivlg II region, and the radiative transfer of the lines 
is very different. Computing the Fe II and JVIg II line 
strengths for a range of hydrogen densities and ionization 
parameters and for a cosmic abundance of elements, they 
predict values between 1.5 and 4 for the line ratio and 
attribute higher ratios to an overabundance of iron, with 
respect to magnesium, but probably also with respect to 
hydrogen. These computations are, howe ver, based on a 
limited 70-level model of the iron atom. iBaldwin et al.l 
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TABLE 2 

Fitted spectral properties of the FeII and Mgii lines and FeII/Mgii ratio 



QSO 




a'' FfoIi" 


^Mgll'l -pMgll'^ 
10-15 


EWcMgii*^ FWHMMgii Fcll/Mgll 






10-15 






[erg cm— ^ s-l] 


[erg cm— ^ s-i] 


[A] [kms-i] 



J08364 


-0054 


5.82 


-1.81 


7.4dzl.2 


5.810±0.003 


2.15±0.17 


24±3 


3600±300 


3.4±0.6 


J0005- 


-0006 


5.85 


-2.05 


2.3±0.6 


5.850±0.003 


0.48±0.06 


29±5 


2100±300 


4.7±1.5 


J141H 


-1217 


5.95 


-1.96 


2.7±0.7 


5.904±0.002 


1.13±0.07 


28±3 


2400±150 


2.4±0.6 


J13064 


-0356 


5.99 


-1.94 


2.4±0.9 


6.016±0.002 


1.14±0.04 


35±2 


4500±160 


2.1±0.8 


J10304 


-0524 


6.28 


-1.94 


1.2±0.3 


6.308±0.001 


0.55±0.01 


25±1 


3600±100 


2.1±0.5 


^ The re 


dshift 


as published in 


the diseove 


y paper, based 


on the observed 


wavelength of the Lya line.l* 


Slope of 



the power law continuum: Fx 
2850 < Ao < 3070 A).'^ 



oc A"."^ Flux derived from the Fell template normalization (derived over the range 



The redshift derived from the Mgll centroid.' 
frame equivalent width of Mgll line in A. 
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'e 



1450 



Rest-frame wavelength (Angstrom) 
1500 650 1600 



1650 




1.10 
1600 



1650 



3.0 



0.0 



z J1411+1217 A 









n- 








1.00 




1.05 




1.10 












1400 




1450 


1500 


1550 




1600 




'< 
B 

u 


3.0 
2.5 




J1306+0356 
















2.0 




















aa 


1.5 
1.0 




















b 


0.5 






















0.0 




























1.00 




1.05 




1.10 








1500 




1600 




1700 


1800 




1900 





2.5 ; 
2.0 ; 



J1030+0524 : 



0.0 ; 



1.10 



1.15 



1.20 1.25 1.30 

Observed-frame wavelength (micron) 



1.35 



1.40 



Fig. 5. — Z— and J-band spectra of the five z ~ 6 QSOs, covering the wavelength region that includes CIV, smoothed over five pixels 
(or 0.0014 /xm). Overplotted is the simultaneous fit of a power law and Lorentzian curve (solid [red] line). The absorption feature within 
the Civ line in J1411+1217 was fitted simultaneously with a Voigt function. The unabsorbed Lorentz profile for J1411+1217 is shown 
with a dashed line. The dashed [blue] line shows the noise level per pixel (smoothing is not taken into account here). Overplotted on the 
spectrum of J1030-f0524 in grey [green] is the telluric transmission measured, showing the problematic telluric absorption at 1.136 /xm. 
[See the electronic edition of the Journal for a color version of this figure.] 
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TABLE 3 

Fitted spectral properties of the C IV line 



QSO 


2=^ 


^Civ*^ F'^ 
10-15 

[ergcm"^ s~^] 


EWo'* 
[A] 


FWHM 

[km s-i] 


J0005-0006 


5.85 


5.848±0.001 


l.SitO.l 


40±2 


2900±140 


J1411+1217 


5.95 


5.911±0.001 


S.OitO.l^ 


70±3 


3060± 80 


J1306+0356 


5.99 


5.997±0.001 


3.9±0.2 


40dz2 


5940±200 


J1030+0524 


6.28 


6.262±0.003 


2.4±0.2 


32dz3 


5030±380 



^ The rcdshift as published in the discovery paper. ^ The rcdshift de- 
rived from the centroid of the CIV Une."^ Flux derived from the CIV 
line fit.*^ Rest frame equivalent width in A.^ Unabsorbed value, actual 
observed value is smaller. 



(|2004f l construct a 371-level Fe+ model, using energy bal- 
ance to obtain a self-consistent temperature and ioniza- 
tion structure. They apply this model to a large range of 
parameters in order to reproduce the observed Fe II emis- 
sion properties and the relative strengths of the strong 
emission line of other elements (such as Lya and Mg II) . 
Only two scenarios for observed quasar spectra are al- 
lowed by the results of this modeling: there is either 
significant microturbulence (fturb > 100 kms~^), consis- 
tent with currently popular models where the emission 
lines are formed in a wind flowing off a rotating accre- 
tion disk, or the Fe II lines are formed in gas with other 
properties than the gas where the other emission lines are 
formed, i.e., in spatially separated regions. The bottom 
line is that, although the relative stren gth of Fe II does 
depend somewhat (but not hnearly, Verner et al.lll999h 
on the iron abundance, it also depends sensitively on 
other parameters, preventing strong conclusions on the 
Fe/Mg abundance ratio from the Fe II/Mg II line ratio, 
including the comparison with the presumed solar value 
(see Fig. [6]). 

Although our Fe II/Mg II l ine ratios are consis- 
tent with the va l ues d erived by llwamuro et al.l ()2004l ) 
and iBarth et al.l (|2003[ ). t hey are systematicall y lower 
than those measured by iMaiolino et al.l (120031 ). For 
J1030-I-0524 ou r ratio is consistent with that found by 
iFreudling et"all (|2003f) . while for J0836-t-0054 it is lower. 
Note that, using the o r iginal Fe II template provided by 
IVestergaard fc Wilkes! ()2001l ) would result in even lower 
Fe II/Mg II ratios for our sample (as we would have 
overestimated the contribution of the Mg II line emis- 
sion). Within individual studies, the emission line ratios 
for different z ~ 6 quasa rs also span a large range of 
val ues (e.g., 1.0 to 9.5 for llwamuro et al] and 2.1 to 9.5 
for IFreudling et all ) . indicating that there is significant 
scatter within the population of z ~ 6 quasars. There 
are several possible causes for the observed discrepancies 
between different authors, including the iron template 
employed, the wavelength limits over which the template 
is integrated, the wavelength range over which the tem- 
plate is actually fit, and the length and signal-to-noise 
ratio of the spectra used. To derive clear conclusions on 
the evolution of the Fe II/Mg II line ratio, a better com- 
parison study should be made where all spectra available 
in the literature are fit by the same template using the 
same method, but this work is beyond the scope of this 
paper. 

We have found no strong correlations between 
Fe II/Mg II line ratio, continuum power law slope a 
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Fig. 6. — Fe II/Mg 11 line ratio plotted versus redshift of our 
sample of quasars (encircled stars) compared with quasars at red- 
shifts < z < 6. 4 from the literature: Dietrich et al. (2003, 
diamonds, including 'Dietrich et al. 2002), Freudling et al. (2003, 
squares), Iwamuro et al. (2004, triangles), Thompson et al. (1992, 
circles), Barth et al. (2003, upside down triangle), Maiolin o et al] 
((2003, ellipses). Filled symbols represent sample averages and their 
error bars represent the standard deviation of the averages of the 
measured values. The range of values obtained by modelling the 
Fell and Mgll line strengths of gas with cosmic abundance and a 
range of ionizing continua (2.75±1.25, see Wills et al. 1985) is indi- 
cated by a grey band. See text for details regarding the difference 
in the determination of this ratio employed by various authors. 



and luminosity within our sample. The scatter in the 
Fe II/Mg II line ratio seems to be small in our sample 
and independent of luminosity, although the lowest lu- 
minosity object is also the iron richest. However, this 
object, J0005— 0006, seems to be peculiar, as it h as nar- 
rower UV emission lines (1500 — 2500 kms-^. Fa n et al.l 
12004) than other SDSS QSOs and may not be represen- 
tative of the general z 6 QSO population. 

We here summarize the (in) consistencies in measure- 
ments by different authors for individual objects in our 
sample. 

Notes on individual objects 

J0836 +0054 This QSO , at z = 5.82 was observed in th e 
NIR bv IFreudling et all (|2003f) and IStern etld] (|2003f) . 
IFreudling et al.l (|2003l ) find a (scaled ) Fe II/Mg II line 
ratio of 9.5±2.3, while the spectrum of lStern et al.l ()2003l ) 
does not include the ii'-band. We find a ratio of 3.4zfc0. 6, 
significantly lower than that found bv IFreudling et all 

J1306+03 56 This QSO, a t z = 5.99 was observed in 
the NIR bv IMaiolino et al.l ()2003l ). who found a ratio of 
9.0±2.3, much higher than the value of 2.1±0.8 found by 
us. 

J1030+0524 The most distant QSO of our sample, at 



6.28 w as observed in t he NIR by 

20031) and iPentericci eFall (|2002| ), 



Maiolin o et al 
Maiolino et al 



mm . again, find a much higher ratio of 8.7±2.5 than 
ours of 2.1±0.5, which, ho wever, is consistent wit h the 
value of 2.1±1.1 derived bv IFreudling et all (|2003f l from 
a NICMOS spectrum sampling rest-frame wavelengths 
between 1600 and 2700 A. The ratio derived by us is 
also consistent with the ratio of 0.99lg gg reported by 
llwamuro et"all (|2004f l. 
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The enormous energy output of QSOs is most probably 
powered by a central black hole which is surrounded by 
an accretion disc. The mass of this black hole is a fun- 
damental parameter which, at least in nearby galaxies, 
is strongly correlated with the bulge velocity dispersion 
of the galaxy host (the so called Mbh — Cbuige relation, 
where (Tbuige is the velocity di s persion of the host galaxy's 
bulge e.g.. iGebhardt et alJ I2OOOI: iFerrarese fc MerrittI 



[25oa iTremaine et al.l 120021: 1 Woo etal.1 120061 ). This re 
lation suggests that the evolution of the black hole and 
the build-up of the bulge in its host galaxy are in- 
timately related. The central black hole mass of our 
Galaxy was measured by very high spatial resolution 
imagi ng and subsequen t modelling of the stellar dynam- 
ics ( Schodel et al.ll2003[ ). but this method cannot be ap- 
plied to extragalactic black holes as the required spatial 
resolution cannot be attained. The most direct measure- 
ments of central black holes mass i n AGN come from re- 
verbaration mapp ing studies (e.g.. iPeterson et all 120041 : 
iKaspi et al.ll2007l ). but these time intensive observations 
are not easily feasible for very distant QSOs. Here, we 
derive black hole masses using three different approaches, 
as described below. 

4.2.1. Eddington mass 

The simplest method to estimate the black hole mass 
is to use the bolometric luminosity of the QSO and to as- 
sume that the outward radiation pressure gradient acting 
on the gas is just enough to counterbalance the gravita- 
tional attraction of the black hole. This gives an esti- 
mate of the black hole mass, called the Ed dington mass 
(|Eddingtonlll922l : iMargon fc Ostrikeij|1973[ ) 



LEdd = 1.3xl03«f^)ergs-i 



M, 







(2) 



i.e., the black hole mass -A/bh can be derived assuming 

4.2.2. Black hole mass derived from Mg II 

Under the assumption that the dynamics of the broad 
line region are dominated by the gravity of the central 
black hole, an estimate of the black hole mass is given by 
Mbh ~ G~^RblrV'^, where Mbh is the black hole mass, 
^BLR the radius of the BLR and V the velocity of the gas, 
which can be estimated fro m the width of emission lines 
in the QSO's spectrum (see iMcLure fc Dunlop 2004 and 
references therein for a discussion of this assumption and 
of the emission lines which can be used). The radius of 
the BLR is strongly correlated with th e AGN monochro- 
matic continuum luminosity at 3000 A (jMcLure fc Jarvii 
l2002f) . allowing black hole mass estimates for high red- 
shift quasars to be made from a single spectrum covering 
the blended Mg IIA 2795,2803 line doublet and redward 
c ontinuum. 

IMcLure fc ( pOOi ) use 17 AGN (aU at z < 0.7) 

with luminosities comparable to those of known high red- 
shift QSOs (ALa > 10^'^ W) for which reverberation 
mapping es timates of t he BH mass are a vailable (from 
IWandel e~ al. 1999, and lKaspi et al.l[2000[ l. to determine 
the following relation between black hole mass, contin- 
uum luminosity at Aq — 3000 A and Mg II line width: 



BH 



M^ 



= 3.2 



[Aisooo "I 
1037 w J 



f FWHM [Mgll] 



kms~i 



(3) 



where L3000 is the luminosity per wavelength bin at rest- 
frame Ao = 3000 A and FWHM is the Full Width at Hal f 
Maximum of the Mg II line. iMcLure fc JarvisI ()2002D 
state that an earlier version of this black hole mass esti- 
mator can reproduce the reverberation black hole mass 
to within a factor 2.5 (Icr). 

4.2.3. Black hole mass derived from C IV 

IVestergaard fc Peterson! ()2006D use 32 AGN with re- 
liable reverberation mapping estimates calculated by 
IPeterson et al.l ()2004D to derive the following empirical 
relation between black hole mass, continuum luminosity 
at Ao = 1350 A, and C IV hne width: 



M, 



BH 



M^ 



= 4.6 



[ ■^-^1350 ^'^'^ ( 
1037 W J [ 



FWHM [CIV] 1 



kms 



J 



(4) 



IVestergaard fc PetersonI estimate that this mass estima- 
tor is accurate within a factor 3.6 — 4.6, but they also 
note that the reverberation-based masses themselves are 
uncertain typically by a factor ~ 2.9. 

It should be n oted that some authors (e.g., 
iBaskin fc Laoil 120051 ) have argued that the C IV orig- 
inates, at least in part, from outflowing gas or that 
the shift and asymmetry commonly seen in C IV sug- 
gest that non-gravitational effects, such as obscuration 
and radiation pressure, may affect the line profile (e.g., 
iRichards ct al. 2002^) and its FWHM would therefore 
not be a good measure of the gas dynam ics in the BLR, 
although IVestergaard fc PetersonI (|2006f l claim that the 
typical uncertainties resulting from these are smaller 
than those quoted above for the mass estimator itself. 

4.2.4. Derived black hole masses 

Using these three approaches, we estimate black hole 
masses for our sample and list them in Table 21 The 
bol ometric lum inosities for these quasars were obtained 
by iJiang et al.l (2006.) from SED template fits to com- 
bined observations at wavelengths including the radio, 
mm/submm, UV/optical, NIR, MIR and X-rays. For 
the computation of ^1359, we had to extrapolate the con- 
tinuum measured to the wavelength required using the 
fitted power law, but given the fact that the continuum 
is well characterized in the Z- and J-band, this should 
not introduce large errors. The black hole masses for the 
sample at hand range from 0.3 to 5.2 x 10^ Mq, includ- 
ing the lowest black hole masses ever observed for z 6 
QSOs. In general, the black hole masses derived using 
the three different techniques agree within a factor of a 
few, where the mass derived from the C IV line is largest 
and the Eddington mass in between. We note that, even 
in the absence of any measurement errors, the individ- 
ual BH mass estimates have an uncertainty of a factor 
^ 3 from the scat t er in the c orrelations measured by 
IMcLure fc DunlopI (|2004( l and IVestergaard fc PetersonI 
(2001). 

If the quasars were strongly lensed, the mass de- 
rived by the Eddington luminosity would be larger by 
y/g, where q is the lensing factor, as Msdd ~ q and 
-^'^MgH/civ roughly scales with ^/q). Assuming that 
qu asars emit isotropically (whi ch may not be the case, 
see lHennawi fc Prochaska 2007*). the good agreement be- 
tween the different black hole mass estimates implies that 
the quasars are likely not strongly lensed, consistent with 
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TABLE 4 

Estimated black hole masses 



TABLE 5 
Redshift measurements 



QSO 




A/BH(MgIl) 


AfBH{ClV) 


AfBH(Edd) 






[109 Mq] 


[109 Mq] 


[109 Mq] 


J0836+0054 


5.82 


2.7±0.6 




4.0 


J0005-0006 


5.85 


0.3±0.1 


0.7±0.1 


0.7 


J1411+1217 


5.95 


l.litO.l 


1.2±0.2 


1.2 


J1306+0356 


5.99 


2.4±0.4 


5.2±0.7 


1.9 


J1030+0524 


6.28 


1.4±0.2 


3.2±0.6 


1.8 



^ The redshift as pubhshed in the discovery paper. 

the conclusions derived from an analysis of HST images 
of these objects, non e of the QSOs showing eviden ce for 
lensed counterparts (iRichards et al.l |2004 l2006a[ ). For 
J1030+0524, Haiman & Cen (200^ predict a minimum 
BH mass of ^ IQ'^Mq based on the flux distribution of 
the Lya emission alone and assuming it is not gravi- 
tationally lensed. The BH mass derived by us for this 
quasar is in the r ange 1.4 — 3.2 Wi^Mp^ , consistent with 
the prediction bv iHaiman fc CenI ([2002) . 

Black hole masses derived from the observed-frame 
NIR emission hues for the QSOs in our sample have 
not been reported in the literature, although for 
some of the sources t he C IV line h as been ob- 
ser ved (J0836+0054 bv |Stern et all [2001 J1306+0356 
by iMaiolino e t al. '2004a' and 'Pc ntericci et al.l l200l 
J1030+0524 by Maiohno ct al. 2004ai ). We can compare, 
however, with the BH mass measured for the z = 6.4 
QSO J114 8+5251, which was measured by iBarth et al.l 
(|2003f ) and lWillott et all (|2003f l. employing relations be- 
tween the Mg II and C IV line widths and black hole 
masses very similar to those employed by us. These au- 
thors find a BH mass of (2 - 6) x 10^ Mq and 3 x 10^ 
Mq, respectively. The BH masses measured by us are 
in general smaller, showing that we sample the BH mass 
function to lower masses than previous studies, probing 
not only the most massive and luminous objects but also 
the presumably more common less massive quasars. 

4.3. Redshifts of QSOs in our sample 

For several reasons, it is useful to know the systemic 
redshift of a QSO with high precision, e.g., for follow-up 
observations of molecular lines with narrow bandpasses. 
The systemic redshift is difficult to derive from the cen- 
troid of the observed Lya line, as it is very suscepti- 
ble to absorption by neutral hydrogen on its blue side, 
can be blended with the N V line or because it is com- 
pletely absent. Narrow emission lines due to star forma- 
tion are more suitable redshift indicators. For example, 
the [O III]A5007 line is presumed to be at nearly the 
systematic redshift of th e host galaxy (see Sec. 3.0 of 
IVanden Berk et all 1200 iL for references). If narrow lines 
cannot be observed, broad QSO emission lines other than 
Lya can give important information about the systemic 
redshift of the galaxy obser ved. Using the red s hift o f 
[O III]A5007 as a zeropoint, IVanden Berk et all ()200lD 
find an average redshift offset of 143±91 and 161±10 km 
s^^ for Lya and Mg II and a blueshift offset of 563±27 
km for C IV (where the errors indicate uncertainties 
i n the mean valu e , the i ntrinsic scatter is much larger). 

IRichards et al.l (j2002b) study these offsets in more de- 
tail and find a redshift offset of 97 km s^^ for Mg II and 
a blueshift offset of 727 km s^^ for C IV, with dispersions 
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2^ 


zciv 


ZMgll 


, b 
Zsys 


J0836+0054 


5.82 




5.810±0.003 


5.808±0.007 


J0005-0006 


5.85 


5.848±0.001 


5.850±0.003 


5.848±0.007 


J1411+1217 


5.95 


5.911±0.001 


5.904±0.002 


5.902±0.007 


J1306+0356 


5.99 


5.997±0.001 


6.016±0.002 


6.014±0.007 


J1030+0524 


6.28 


6.262±0.003 


6.308±0.001 


6.306±0.006 


Sv'^ (km s"-'-) 




-193 


97 




5z<^ 




0.005 


-0.002 





^ The redshift as pubhshed in the diseovery paper. ^ Systemic redshift 
estimate based on corrected Mgll redshift.'^ Velocity difference between 
the listed line and the [OIII]A5007 centroid ffrom lVanden Berlc et al.l 
[2OO1I for Lya and N V, from IRichards eTall I2002bl for C IV and 
Mgll).'^ Correction in redshift for z — 6.0. 



of ±269 and ±511 km s"^, respectively. Their analysis 
shows that the apparent shift of the C IV line peak is 
not so much a shift as it is a lack of flux in the red wing. 
They also confirm the observed anticorrelation between 
the shift of the C I V line peak an d the rest equivalent 
width of this line. IRichards et al.l present this correla- 
tion as follows: they first sorted their list of 794 quasars 
according to their C IV— Mg II redshift differences and 
then divided this list into four roughly equal bins with 
~200 quasars each. The four bins, having C IV EWqS of 
30.3, 25.8, 22.3, 18.0 A, have average blue shifts of 193, 
651, 1039, 1596 km with respect to the Mg II line 
peak, respectively. As all our measured C IV EWqS are 
larger than 30 A, we assume the C IV line is blue shifted 
by 96 km with respect to the systemic velocity. As 
the shifts of line peaks of emission lines in the spectrum 
of an individual quasar are correlated, one can obtain 
a better redshift correction if the peak wavele ngths of 
several lines are measured (see lShen et al.|[2007l ). but we 
do not apply this method as we do not have enough suit- 
able lines with a high enough signal-to-noise ratio in our 
spectra. 

Table ini presents for each QSO, the observed C IV and 
Mg II redshift and the systemic redshift Zgys, which is 
equal to the Mg II redshift corrected as described above 
(i.e., +97 kms~^). The error on the sy stemic redshift 
takes into account the dispersion found bv lRichards et all 
([2002b) for the offset of Mg II (i.e., ±269 kms"!). Rela- 
tively large inconsistencies remain between the corrected 
redshifts of C IV and Mg II. One cause could be that 
C IV, being a resonant line, can potentially suffer from 
absorption, as is the case for J1411+1217 (see Fig. [5]). 

Notes on individual objects 

J0836+0054 This QSO was obser ved by I Stern et ap 
12003). From the C IIIA1909 line, IStern et al.l ([20031 ) 
derive a corrected'* redshift of z = 5.774 ± 0.003. This 
is quite differ ent from our systemic redshift of Zgys = 
5.808 ± 0.007. iFreudling etall ([2003^ measure an uncor- 
rected C IIIJA1909 and Mg II redshift of z = 5.82 from 
their low resolution NICMOS spectrum. 

J1306+0356 Measure ments to the sa me data presented 
here were published bv i Pentericci et a l. (2002), who find 
an uncorrected C IV redshift of z = 6.00±0.01, consistent 
with our result and with the value derived fom N V by 



Usin g the 224 km s~\ veloci ty difference with [OIII]A5007 



U sin g the 224 km s velocit 
found bv iVanden Berk et al.l 1I2OOII ') 
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6500 7000 7500 _ 8000 8500 

Observed-frame wavelength (Angstrom) 

Fig. 7. — Optical spectra of the quasars in our sample, off- 
set in the vertical direction for clarity by 1.5, 4.5, 9.5, and 11.5 
xlO"^'^ ergs~^ cm~^ A~^. The size of the Stromgren sphere (-Rs) 
is indicated by the line which connects the measurements of the 
onset of the Gunn-Peterson (GP) trough (zgp) the systemic 
redshift (zsys) derived from the Mgll line. Error bars for both the 
GP trough onset and the systemic redshift have been indicated, 
although the uncertainties on the latter are too small for the error 
bars to be discernable for most quasars. 



iBecker et all ()2001[ ). 

J1030+0524 iFreudling et"al] (|2003[ ) measure an un- 
corrected C III]A1909 and Mg IIA2800 redshift of 
z = 6.28 from t heir l ow resolution NICMOS spectrum. 
Ilwamuro et al.l ()2004l ) find an uncorrected Mg II red- 
shift of z = 6.3 1 1 (and do not report the uncertainty). 
iPentericci et al.l ()2002[ ) derive uncorrected redshifts of 
z = 6.29 ± 0.01 and 6.27 ± 0.02 from N V and C IV, 
respectively, while we find a consistent redshift of z = 
6.262 ± 0.003 (uncorrected) from C IV. 

4.4. The size of the ionized sphere around the QSO 

At z ~ 6, most of a quasar's emission below the Lya 
line is absorbed by neutral hydrogen in the IGM. How- 
ever, the quasar is able to ionize a (Stromgren) sphere 
around it, from which it transmits radiation at wave- 
lengths below that of the Lya line. The red edge of the 
GP trough marks the radius of this sphere. Therefore, 
the radius of the Stromgren sphere can be derived from 
the difference between the redshift of this red edge and 
the systemic redshift. However, if the GP effect is not 
complete, the precise redshift zqp of the red edge may 
be difficult to determine and may depend on the spectral 
resolution of the observations. We have tried to deter- 
mine zgp conservatively by eye, determining where the 
Lya emission from the QSO may be fully absorbed, but 
including in its error the range of redshifts where this 
plausibly could happen. This results in relatively large 
error bars for some QSOs in our sample, as shown in 
Fig. [T] As the uncertainty in zgp is much larger than 
the uncertainty in the systemic redshift, we disregard 
the latter in the following computation. 



The Stromgren sphere radius {Rs) is e qual to the phys- 
ical distance between Zgys and zgp (see iHog jll999| l: 



i?s = 



(1 + ^) 



\^l^{l + zf + n„)-^'^dz, (5) 



where Dh is the Hubble distance and and J^a the 
ratios of the mass and vacuum energy density to the 
critical density, respectively. For the cosmological pa- 
rameters used here and a redshift z = 6.0, we obtain 
i?s = 60 (zgys — Zqp) Mpc. The derived values of Rs for 
our sample are listed in Table [6l 

4.5. QSO activity time scales 

The size of H II regions around a luminous quasar, 
prior to complete reionization, depends on the ionizing 
photon emission rate (TVph), the time span the ionization 
is sustained (tq) and the density of neutral hydrogen 
(nna^Hi) in which the sphere expands 



Rs 



3NphtQ 



1/3 



(6) 



where xhi is the neutral hydrogen fraction. This formula 
results from matching the total number of ionizing pho- 
tons emitted to the number of hydrogen atoms within a 
sphere, and it i gnores both recom bina tions and the Hub- 
ble expansion (jWhite et al.ll2003l . see lMaselh et al.ll2007l 
for a discussion of the validity of the assumptions made 
here). If the quasar is located in a region overdense by 
a fac tor D, the radius is reduced by D~'^/^ (jWhite et al.l 
l2003fl . Furthermore, from the analysis of mock spectra 
along lines of sight through a simulated QSO environ- 
ment, iMaselli et al.l ()2007) find that the H II region size 
derived from quasar spectra is on average 30% smaller 
than the physical one. Note that issues of geometry, H I 
clumping factor, and pre-QSO reionization by, e.g., star 
forming galaxies in the QSO vicinity, further complicate 
a conversion from observed Stromgren sphere radius to 
quasar activity time. As we are interested only in or- 
der of magnitude effects, we do not take these issues into 
account. 

Assuming that the quasar is located in a completely 
neutral IGM (i.e., xm ~ 1) which contains nearly all the 
baryon s (i.e., rtn is equal to the r nean cosmic density at 
z = 6'l. lCen fc Hah^ ((2OO0I) and lHaiman fc CenI (poM ) 
use the the radius of the Stromgren sphere to constrain 
the minimum activity time for the quasar 



to = 



A^t 



ph 



1057 s- 



Rs 



7.4 Mpc 



7.0 



lOVr- (7) 



Note that this time is very short, and of th e order of the 
light c rossing time of the Stromgren sphere. iWhite et al.l 
(|2003f ) explain that the speed of light drops out of the 
equations describing the early evolution of a Stromgren 
sphere, and that, for an observer, there is an initial period 
of superluminal expansion. 

The ionization rate depends on the SED below (blue- 
ward of) the Lya line, which we estimated from the SED 
above 1216 A. The composite quasar SED of ' Elvis et all 
(I994) with median bolometric luminosity lO^'' erg s~^ 

has a corresponding median Aph — 5.5 x 10^'^ photons 
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TABLE 6 

Stromgren sphere radii and quasar activity time scales 



QSO 






^GP 


Zsys 


Rs 
[Mpc] 


.^^ph.EIv 

[1057 s-1] 


*Q.EIv 

[10^ yr] 


*Q.TcI 

[10^ yr] 


J0836+0054 


5.82 


-27.83 


5.721±0.033 


5.808 


5.7±2.2 


4.70 


0.6 


0.1 


J0005-0006 


5.85 


-26.42 


5.771±0.044 


5.848 


4.9±2.8 


0.47 


3.9 


0.8 


J1411+1217 


5.95 


-26.70 


5.823±0.021 


5.902 


5.0±1.4 


1.07 


1.8 


0.4 


J1306+0356 


5.99 


-27.14 


5.916±0.059 


6.014 


5.9±3.6 


1.40 


2.4 


0.5 


J1030+0524 


6.28 


-27.10 


6.217±0.020 


6.306 


4.9±1.1 


1.35 


1.4 


0.3 



^ The rcdshift as published in the discovery paper. ^ AB magnitude at rest-frame wavelength 1450 A 
from the discovery papers, corrected for the cosmology used here.*^ Rcdshift from the centroids of the 
CO lines measured by [Walter et al.ll2003l and lBertoldi et al.ll2003bt 



s~ . We have used this relation to derive the ioniza- 
tion rates for our sample QSOs (A'^ph^Eiv, see Tab. [5]). 
From these ionization rates and Eq. [7] we derive mini- 
mum activity times for the quasars in the range: ^q.eiv ~ 
0.6-9x10"^ yr. However, as noted by|miteet al. (2003), 
who h ave used the quasar te mplates by [Telfer et al.1 
()2002f ). the ionization rates bv lElvis et al.l may be too 
low by a factor of five. Consequently, the quasar activity 
times derived from the Stromgren spheres would go down 
by the same factor: <Q,Tei ^ 0.1 — 2 x 10^ y r (see last col- 
umn of Tab. O- For a partly {xm =0.1. IWvithe et aTl 
[20051) or even mostly (xhi = 0.01, iMaselh et alj 1200*^ 
ionized IGM at z ~ 6, the derived activity times are 
smaller by a factor 10 or even 100. 

We can now compare these ages to the e-folding time 
scale for t he central accr eting black hole (the Salpeter 
time scale. ISalpeter|[l96^ 

tacc = 4x 10^(^)7,-1 yr (8) 

where e is the radiative efficienc y and rj is the ratio of 
quasar to Eddington luminosity (Haiman fc Loebll200lL 
see also the discussion in .White ct al. 2003). Assuming 
that the quasar is radiating at the Eddington limit (77 ~ 
1) and an efficiency of e ~ 0.1, this simple calculation 
gives values that are consistent, to first order, with our 
derived ages of the Stromgren spheres in a neutral IGM. 
This implies that the mass of the central black hole in 
the quasars observed grows only by one e-folding or less 
during one quasar activity cycle. 

4.6. Absorption features 

No Broad Absorption Lines (BALs) are found in the 
spectra, although the C IV emission line of J1411-I-1217 
shows associated C IV absorption. The blue shifted 
BALs in QSOs are probably caused by an outflow- 
ing wind. At low and intermediate rcdshift, BAL 
QSOs account for ~15% of the whole quasjir popula- 
tion (iReichard et al.l 120031 : iHewett fc Foltzl l2003f l. but 
llVraiolino et all (|2004a[ ) ffiTd four BAL QSOs among the 
eight QSOs at 4.9 < z < 6.4 that they observe, sug- 
gesting that the most distant quasars are surrounded by 
a much larger amount of dense gas than lower rcdshift 
quasars. Adding the three QSOs without C IV BALs of 
our sam ple (J1030-I-052 4 was already found to be non- 
BAL bv iMaiohno et al.1 ). the fraction of BAL QSOs at 
z > 4.9 (5.7) goes down to 36% (25%). A larger sam- 
ple of z ~ 6 QSOs with C IV spectroscopy is needed to 
resolve this issue. 

While not immediately the science focus of this pa- 
per, we have also searched for narrow absorption lines of 



intervening systems. The spectral resolution of ~ 20 A 
in our sp ectra does not allow to identify the C IV dou- 
blet (see ISimcod 120061: iRvan- Weber et al.ll2006D . but it 
is good enough to resolve the redshifted lines of the 
Mg IIA2796,2803 system. Within the continuum emis- 
sion of the QSO J1306-I-0356, we find two systems of 
double absorption lines: at 9797.7, 9827.4A and 9880.7, 
9906. 3 A with observed EWs of 6.2, 4.0, and 11.6, 10.3 A, 
respectively. In addition, a third system of double ab- 
sorption may be present at the red edge of the J-band 
spectrum: at 11340.6, 11401. 2 A with observed EWs of 
13, 11 A, respectively (see Fig. |5]). These absorption 
features would be consistent with Mg II absorption at 
z = 2.504, 2.533 and 3.060 (although the latter only 
marginally), but a systematic search with larger wave- 
length coverage by Jiang et al. (2007, AJ, in press) 
shows that the correct identifications for the first dou- 
ble absorption feature is Al II 1670 at z = 4.864 and 
4.882, while the third double absorption feature may be 
Fe III 1926 from the same intervening systems. We note, 
however, that the latter lines are at the edge of the ob- 
served spectrum where the noise is higher and the wave- 
length calibration is not very accurate. The second dou- 
ble absorption feature is indeed confirmed to be Mg II 
at z = 2.533 from other absorption lines identified. The 
two absorbers at z = 4.9 a re the highest rcdshift Mg II 
absorbers known (see, e.g.. iKobavashi et al.l 120051 ). See 
Jiang et al. for more details and a discussion on the 
absorber host galaxies. 

4.7. Continuum slope 

The continuum slopes are measured over small wave- 
length intervals with little pure continuum emission, 
especially for the A'-band spectra. Despite result- 
ing uncertainties, the measured slopes lie in the range 
—2.06 < a < —1 - 83, st eeper than that measured by 
IVanden Berk et all (|200lL -1.54) for the SDSS composite 
QSO spectrum at a median 2 = 1.3 and than that mea- 
sur ed from NIR pho t ometr y of 45 QSOs at 3.6 < z < 5.3 
bv JPentericci et al.l (|2003L -1.43±0.33). iMaiolino etlHI 
(|2004a| ) find such steep continuum slopes for some of 
their QSOs at 4.95 < z < 6.40, and in particular for 
J1048-I-4637 at z = 6.19, where they measure a slope 
of g = —2.1 at rest frame wavelength above 1700 A 
(jMaiolino et al.ll2004bD . Jiang et al. (2007, AJ, in press) 
observe a z ~ 6 QSO sample with larger wavelength cov- 
erage and are therefore able to constrain the continuum 
slopes better, obtaining results similar to ours. 

5. SUMMARY AND CONCLUSIONS 
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with the conclusions derived by previous studies (see also 
Barth et all 120031 : iMaiolino et all [200l iFreudling etal\ 



1.00 1.05 1.10 

Observed-frame wavelength (micron) 

Fig. 8. — Z-band spectrum of J1306+0356, smoothed over three 
pixels, with AlII absorbtion from systems at 2 = 4.864 and 4.882 
and Mgll absorbtion at z = 2.533 indicated. The grey [green] line 
shows the noise per pixel. [See the electronic edition of the Journal 
for a color version of this figure.] 

We present sensitive near-infrared spectroscopic ob- 
servations of a sample of five z ~ 6 quasars, compris- 
ing all published quasars at z > 5.8 accessible from the 
VLT. Our ISAAC spectra cover the C IV, IVIg II and 
Fe II lines, which are powerful probes of the chemical 
enrichment and the black hole masses in these objects. 
As the Mg II and Fe II lines are blended (and are lo- 
cated on top of the Balmer continuum), we have used 
Fe templa tes derived from an SPS S composite quasar 
spectrum (IVanden ^ Berk et"al?'2001') and from the spec- 
trum of I Zw 1 (jVestcrgaard & Wilkes 2001) and simple 
models representing the continuum and Balmer pseudo- 
continuum to decompose the spectra. We derive an av- 
erage Fe Il/lVIg II (which is a proxy for the Fe/a) ra- 
tio of 2.7±0.8 for our sample, which is similar to the 
values measured for quasar broad line regions at lower 
redshifts. We note that the values measured by us are 
lower by a factor of a few than some previous mea- 
surements for ind ividual quasars at z > 5.8 reported 
in the literature fJVIaiolino et al.l l2003t IFreudling et al.l 
[2003: Iwamuro et al. 2004) but are consis t ent with others 
(iBarth et al.ll2003t iFreudling et al.ll2003l : llwamuro et all 
|2004| ). We attribute these differences to different choices 
in the Fe templates used, the corresponding spectral de- 
composition and the different quality of the spectra. The 
lack of evolution in the observed Fe II /Mg II ratio and 
the mere presence of iron in the z 6 systems demon- 
strate that the quasars in our sample must have under- 
gone a major episode of iron enrichment in the short time 
span between the Big Bang and the emission of radia- 
tion observed by us, i.e., in less than one Gyr, consistent 



We have derived central black hole masses using three 
different methods: A) using the Eddington luminosity, 

B) using the empirical equation relating the linewidth 
of IVIg II (and the nearb y continuum emiss i on) t o the 
mass of the black hole (|McLure fc Dunlod I2004D and 

C) using the corresponding m ethod for the CIV line 
(jVestergaa rd fc Petersonll2006f ). The derived masses us- 
ing the JVIg II and C IV lines and the Eddington luminos- 
ity agree well, and we derive central black hole masses of 
0.3 — 5.2 X 10^ Mq. Our derived black hole mass range in- 
cludes the lowest black hole mass ever measured at z ~ 6 
and indicates that we have started to probe the more typ- 
ical quasar population (with lower masses/luminosities) 
at this extreme redshift. Clearly, future, more sensitive, 
observations of an extended sample towards lower lumi- 
nosities are needed to constrain the properties of the 
fainter quasar population further. Assuming that the 
quasars indeed radiate at their Eddington luminosity, 
the agreement in derived black hole masses using the 
different methods implies that the quasars are likely not 
strongly lensed and emit roughly isotropically. 

We use the difference between the redshift derived 
from the Mg II line, which serves as the best available 
proxy for the systemic redshift of the quasar, and the 
redshift for the onset of the CP trough in the quasar 
spectra to derive the extent of the ionized Stromgren 
spheres around our target quasars. The derived physical 
radii are relatively well constrained to about five IVIpc. 
Using simple ionization models, the central quasar 
activity would need of order 10^ — 10* yr to create these 
cavities, allowing the quasar black hole to grow by one 
e-folding or less during an activity cycle. 
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